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Abstract 

The correct Hamiltonian for an extended Hubbard model with quantum 
group symmetry as introduced by Montorsi and Rasetti is derived for a D-dim- 
ensional lattice. It is shown that the superconducting SUg(2) holds as a 
true quantum symmetry only for D = 1 and that terms of higher order in 
the fermionic operators are needed in addition to phonons. A discussion of 
Quantum symmetries in general is given in a formalism that should be readily 
accessible to non-Hopf algebraists. 
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1 Introduction 



The Hubbard model is the simplest description of itinerant interacting electron 
systems. In this article we will study generalizations of it on a D-dimensional lattice. 
The Hamiltonian of the standard Hubbard model is given by [Q] : 

i7^„, = if(p-^-) + if(;-) (1.1) 

where 

H^-n-loc)^_^ E (1-2) 

<j,j>,(7 

i i,a 

The 1-dimensional model has been solved in [0 . It is well-known that the Hubbard 
model has a (SU(2) x SU(2))/Z2-symmetry This symmetry is the product of 

two separate SU(2) symmetries: a magnetic and a superconductive one. 

Montorsi and Rasetti [^] have introduced a very interesting generalization of the 
Hubbard model by adding phonons. It turns out that the symmetry of the stan- 
dard Hubbard model is sometimes a special case of a more general quantum group 
symmetry. More precisely, while the "magnetic" SU(2)-symmetry is left unchanged, 
Montorsi and Rasetti claimed that the generators of a "superconductive" SUg(2) 
quantum group commute with their extended Hamiltonian. We were able to verify 
this symmetry for an extended Hubbard model on a one-dimensional latticeQ while 
we found unsurmountable obstructions in the higher dimensional case. As we will 
show this is essentially due to ordering problems. Our task in this article is twofold: 
we will address quantum symmetries in general and we will carefully re-examine the 
Hubbard model with phonons, deriving each term on physical grounds to obtain the 
correct Hamiltonian. 

2 Quantum symmetries in quantum mechanics 

The role of symmetries in quantum mechanics cannot be underestimated. Some 
models (harmonic oscillator, hydrogen atom) were in fact first solved relying only 
on symmetries. Symmetries, especially infinite dimensional ones, serve to provide 
the constants of motion that are central to integrable models. 

It is interesting to see what happens when the usual notion of symmetry is relaxed 
and transformations given by a Hopf algebra (Quantum Group) are considered. 
^ The Hamiltonian in pj is given explicitly only in the one-dimensional case. 
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To simplify the discussion we will use a formalism that avoids direct reference to 
Hopf algebraic methods. As given data we take a *-Hopf algebra W, its dual Hopf 
algebra U* and a *-algebra A generated by quantum mechanical operators that act 
on a Hilbert space Ti. The generators of quantum symmetry transformations live 
in U. Here we typically have a one or more parameter deformation of the universal 
enveloping algebra of a Lie algebra in mind. The elements of the dual Hopf algebra 
U* play the role of functions on the Quantum Group. The only difference to the 
classical case is that these functions no longer commute. 

Unitary representation The elements oiU should act on Ti. We need a unitary 
representation p on H that reahzes U in the operator algebra A. Such a represen- 
tation shall be a linear *-preserving map 

p:U^A, p{x):H^H, p{x)^ = p{x*) (2.1) 

that is also an algebra homomorphismQ 

p{xy) = p{x)p{y). (2.2) 

Here is an example: 



Magnetic and Superconductive SU(g)(2) The algebra of SUg(2) is generated 
by X~ = (X+)* and H = H* with deformed commutation relations 

[H, X±] = ±2X±, [X+, X-] = ^^f^; q e R\{0}. (2.3) 

As can be checked by direct computation this algebra has the same representation 
by 2 X 2 matrices as the undeformed SU(2), namely 

-^-{iD- 

From this matrix representation we can find a unitary representation of the algebra 
( P^.3| ) by creation and annihilation operators using the following simple observation: 

Let c|, Ci be fermionic or bosonic creation and annihilation operators 
and rriij, riij numerical matrices with the same (finite) index set as the 
c\ c, then [c^^ ■ m ■ c, ■ n ■ c] = ■ [m ; n] ■ c. 

^Remark: p{x)p(y) — p{z) <^ p{z) — p{xy) <;= z = xy, but not: z — xy" 
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If we take for instance c| G {fej, &|} and Cj G {fef, 6|} and the matrices from ( ^.41 ) we 
find the generators of the "magnetic" SU(2): 

pUX^)=b^% pUX-) = h%, pUH) = {b^%-h%). (2.5) 

Switching ^ &| does not change the algebra of the c^, c (the b\ b are fermionic!) 
but gives another unitary representation — the "superconductive" SU(g)(2): 

p,(X+) = b^%\ ps{X~) = b^b^, ps{H) = (b^^b^ + b^% - 1). (2.6) 

(These expressions hold also for g 7^ 1 because [pm/s{H)]'^ = pm/s{H).) These 
generators implement both (local) SU(2) and (local) SUg(2) for a single lattice site. 
When we deal with generators that act on the whole lattice the "g" reappears and 
consequently (global) SU(2) and (global) SUq(2) do no longer coincide. 
Note: In the sequel we will not write the symbol "p" explicitly; its presence is 
implied by context. 



Transformation of States and Operators The key to a simple description of 
quantum symmetries is the canonical element oiU ®U* sometimes also called the 
"universal T" g 

C = Y,ei®feU®W- (2.7) 

i 

here Cj and /* are (formal) dual linear bases of U and U* respectively. Everything 
else we need to know about C here is that it is invertible and unitary: 

C* =Y,e*® r = C-\ (2.8) 

i 

States IV') G Ti corresponding to a single sit^ transform via multiplication by C: 

|V^)^C|^). (2.9) 

Operators O E A consequently transform by conjugation 

O^COC-^=COC\ (2.10) 

States and operators can have full quantum symmetries, i.e. they can be invariant 
under all of lA. This is the case if respectively: 

^1 /\ I /\ 

COC — 01 (conditions for /u// symmetry) (2-11) 



■^Statements for "single sites" and "multiple sites" of a lattice obviously apply also to a broader 
context of tensor products of states — for instance to single/multi-particle states. 
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When we deal with a lattice, there is a Ci for each of its sites i. Transformations 
of several sites (the whole lattice) i.e. of states \ip^'^^) E 'j-i'^^ and operators O^'^^ E 
^(S>N ^Yso possible. These are performed with products (in the function part) of 
the Ci, 

= C1C2 ...Cn= E e,, e,, . . . (8) e,, ® ■ (2-12) 

11,12, •••i«JV 

SO that 

I^W) ^ c^^)\^i^))^ ^ C(iv)^(iv)(^(jv))-i^ (2.13) 

with (C(^))-i = C^^C^ii . . .Cf ^ Note that the order of the d in C^^^ is important 
because the /' (in the function part of C) are not commutative by assumption for a 
quantum group. 



Full quantum symmetry In the following sections we will be interested in quan- 
tum symmetries of the Hamiltonian. A Hamiltonian h E A has a full "local" sym- 
metry under U (at site i) if 

dhCr^ = /i • 1; (2.14) 
it consequently has a full "global" symmetry under U (on the whole lattice) if 

d^)h(c^^^y^ = C1C2 . ..CNhC^'C^li (2.15) 

In this formalism it is very easy to see that both conditions can also be expressed 
in terms of commutators, namely 



Ci,h 







and 







(2.16) 



respectively. 



Specified transformations Often it is important to describe transformations 
given by specific elements of the Hopf algebra U. So far the transformations were 
unspecific; their result still contained a part in U* i. e. a "function on the quantum 
group" 

e.g. C|^) = Ep(ei)lV') ®r e and similar COC'^eA^W. (2.17) 

i 

A transformation specified by an element k E U is obtained by evaluating these 
function parts on k; this operation will be denoted by "j^". (You may think of it as 
"plugging-in" of the transformation parameters.) The action (denoted by ">") of k, 
on a state l'^) is then given by 

Kt>\ilj) = Cl^l^/j) = k\^Ij) = p{K)\i/j) (2.18) 
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simply because C — being the canonical element — satisfies C|k = ■ = k by 
definition. Similarly 



K>0 



coc-^ 

I K 



(2.19) 
(2.20) 
(2.21) 



The result of contracting the function part of C-^-* = C1C2 ■ ■ ■ Cn with k gives a 
prescription (denoted by /S.^^~^\k) and called the — 1-fold coproduct^) how to 
distribute k over several tensor factors: 



(2.22) 



It is clear that there cannot be one simple rule for all of U — not even in the classical 
case; A(k) = /t®l + l®K for instance holds only for "infinitesimal" n. The added 
difference of the quantum case is that then A(k) will in general be not symmetric. 



Partial quantum symmetry The full quantum symmetries ( p.l4| ) and (|2.15| ) are 
equivalent to 



CihCi 



-1 



h-l^MneU and C^^^h{C 



(7V)n-1 



h ■l^\fK eU 



(2.23) 



respectively. We have seen that these full symmetries could be expressed in terms 
of commutators. As a further illustration of the formalism we will briefiy study 



^The coproduct A did not enter the formalism as additional input here; it rather follows from 
Hopf algebra axioms that 



C1C2 
C1C2C3 



^e,®e,®(/7^')(«;) = ^A(efe)®/'=(^) = A(«), 

i,j k 

(A ® id)A{K) ^{id(g) A) A(k) =: A(2) (k), 



C1C2...C 



N 



The coproducts of a given Hopf algebra are part of the defining relations. Here are the coproducts 
for the generators of the algebra (2.3): 



A{H)=H(g)l + l(g)H, A(X=^)=X^ 



Coproducts of other elements can be computed from this using the fact that A is an algebra 
homomorphism. The other objects that constitute a Hopf algebra are the antipode S and the 
counit e. They enter our formalism via C~^|k = S{k) and 1|k = ^{k). Note that e(K) is a number. 
Let A(k) = ® K(2); then COC-^U - Pie^)OpiSeJ) {fP){K) = p(k(i))Op(5k(2)). This 
action and Hopf-expressions corresponding to equations ( 2.20| -2.21) are e.g. discussed in 



5 



the case where k does not range over all of U but only over a subset V C U. The 
question is: when is 



ChC^^ 
equivalent to 



h-1 



Wk, eV {partial quantum symmetry) 



C,h 



Vk G P 



for an arbitrary Hamiltonian h7 A sufficient condition is easily seen to be 



AC 



VkgP 



(2.24) 
(2.25) 

(2.26) 



for all operators A (g A ^ U*). This can be translated into a condition on the 
coproducts of elements in V: 

A{V)CV®U. (2.27) 

3 A generalized Hubbard model 



Following we will retain the local electron-term (|1.3|) , and add to it the standard 
hamiltonian for the phonons and a phonon-electron interaction term 



Hub 



H. 



(loc) 



+ Hph + 



el—ph- 



(3.1) 



We suppose that the phonons are described by a set of decoupled Einstein oscillators 
with the same frequency u 



(3^2) 



where Pi and jji obey canonical commutation relations as usual. The expression for 
the phonon-electron term is the one given by Hubbard [|l| 

where \E'(r — Ri) is the Wannier electron wave function centered around the ion at 
while fojg., hia are fermionic creation and annihilation operators. (In this context 
the Wannier functions will be approximated by atomic orbitals.) To take account of 
the ion oscillations around their equilibrium positions, the arguments of the Wannier 
functions and of the potential V in the integral must be shifted: 



Rk Rk + Vk, 
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The term obtained from the potential V in has significant contribution only 
foT i = j = I {i.e. neglecting all Ri with I i in V) and results by linear variation 
in a local electron-phonon interaction term ^ 

Hit% = -X-J2in,^ + n,^)y- (3.4) 

i 

and a term that contributes to /i in ( |1.3|) . The non-local electron-phonon interaction 
term is crucial in the approach of Montorsi and Rasetti. We would like to give 
a derivation leading directly to the exponential form, which is necessary for the 
quantum symmetry. (See for the derivation of a linear approximation.) 

We shall retain only the nearest neighbour terms (ij) in the kinetic energy term 
of (|3.3|) ; this assumes negligible overlap between all other atomic orbitals: 

Hir^n'"'^ = E (3.5) 

(iJ) 

with Tij = given by 

T,, = J d^'r m\r- - y,) ^(r - R, - %)• (3.6) 

Assuming that \1' has finite support, it is possible to shift the integration variable 

r — > r — Rj — yj. 

With this substitution Tij becomes a function only of Sij = {Ri + yi) — {Rj + yj): 
T,, = / d^r ^*{f- a,,) (-^] m = T{d.,). (3.7) 



2m 

The atomic orbitals show an asymptotic exponential decay 

^(f) ~ e"^l^^ (3.8) 

and we have hence (approximatly) 

V.,T(a,,) = / d^r *(r-'- a.,) ^ vp(f). (3.9) 

Again, due to the rapid exponential decay of "^{r), we can neglect r in |r — dij\ so 
that 

— * 



V.^,,T(a,,) = -C^T(a,,) (3.10) 

which integrates to: 



l^ij I 



T(a,,) =Toe-'^l'^-l. (3.11) 
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\aij\ = \Ri — Rj + yi — yj\ can be expanded using \yi — yj\ ^ \Ri — Rj\ such that 
finally 

T,, = t exp (^-C ||l|] (yl - y,?j (3.12) 

with a new constant t = Tq exp(— — Rj\). Note that the term 

{Ri - Rj) 



Rij — 



I Ri — Rj I 



always has the same module and that in the one- dimensional case it just amounts to 
a sign. \Ri — Rj\ is the interatomic distance at equilibrium so that it does not depend 
on The complete non-local electron- phonon interaction term in the Hamiltonian 
is 

H^Z'"'^ = ^ E E exp{C4- • iy^ - %•)} b\J^^^ (3.13) 
and the full Hamiltonian of the Hubbard model with phonons is 



Hub 



^E^^T^ii -/^E^i'^ + E + ^^^^Vn - A - E(^iT + ^a)f« 

+ [tJ2 Eexp{C4- ■ im - yj)} + h.cA . (3.14) 

The Hamiltonian considered in ||^ is formally obtained (in the one-dimensional 
case — see remark below) from ( 3.14| ) by a similarity transformation (half of a Lang- 



Firsov transformation ||rT|) on the fermionic operators fo]^ and hia only: 

= U{K)hiM-\K) (3.15) 

a], ^ U{K)h],U-\K) (3.16) 

with a unitary operator 

[/(/3) =exp [i^-Y.Pi^i^ (3.17) 

that depends on a set of constant parameters (3^, k = 1, . . . , D. While this transfor- 
mation does not change the number operators and n^, it results in an exponential 
factor in pi — pj for 

b]A^ = e*'^-®-P^^)a],a,^, (3.18) 
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so that the hopping term is now given by 

H^d-lh^°^^ = ^mi exp{C^ij ■ {Vi - yj)} exp{m ■ {pi - pj)} a]^ai„ + h.c. (3.19) 

{i<j) ^ 

or, combining the exponentials, 

^(non-M ^ Ee"'^'^'^-''exp{C4-(^-%) + ^'?-(p1-p;)} al.a.a + ^-c. (3.20) 

(i<i> 

Remark: Note that while the yi commute with they do not commute with 

the new fermionic creation and annihilation operators a|g., as defined in (|3.15| , 



p.l6| ). The authors of 0, however, assumed commutativity between the fermionic 
operators and coordinates of the ions. In order to be able to connect to their work 
we will formally replace the in i^Hub with new coordinates Xi that do commute 
with the a|^, ai„. (The xi will hence no longer commute with the h\^^ bia-) This will 
of course modify the hamiltonian. The hamiltonian that we will work with in the 
next section is: 

H = H^^""'^ + (3.21) 

with 

= + ^ (A. + \Mu^xi\-\■Y.{n^^+n^l)x,, (3.22) 

i i,a i \ / i 

^{non-ioc) =tY,Y. exp{C^ii " {^i - f,)} exp{m ■ (p, - p,)} a],ai, + h.c. (3.23) 

{i<j) o- 

The relation of this Hamiltonian with the one of the Hubbard model with phonons 
( p. 14 ) will be discussed in section ^. The fact that -ffHub and H are inequivalent can 



for instance be seen by noting that the expression 

fij = texp{(Rij ■ {xi - Xj)}exp{iK ■ {pi - pj)} 



for the hopping amphtude in ( p.23|) does not satisfy the condition Tji = T-j so that 
TijOj^ttia is no longer hermitean. 

4 Superconductive UgSu(2) 

The local superconductive Uq{su{2)) is given by 

p,(X+) = = h\^hl = e-^^-^^^a]^al (4.1) 
p,(X-) = Kt^ = k^k^ = e^*^^^a,^ = (ir;+V (4.2) 
PsiH) = 2Ki'^ = + mi ~ 1 (4.3) 
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These are the generators for transformations of an individual lattice site /, as defined 
in ( p.6| ). They are expressed in terms of the operators bi„, bj^. In order to compute 



the commutation relations with hamiltonian ( 3.21 ) it is necessary to express them 
in terms of operators ai„,aj^ as introduced in ( 3.15| ) and ( ^.16 ). The parameter j3 



appearing in ( |3.17| ), on which the transformation depends, is chosen here to be $/2 



and for the moment it should be regarded as a free parameter which will be deter- 
mined by the commutation relations. We will see later ([4.33|) that consistently with 



the choice (3 = k, made in equations ( |3.15| ) and ( p.l6| ), the commutation relations 



will require $ = 2k. (Notice that n^o- = Oq^CLiu = h\^hifj.) 

To describe the symmetries of the Hubbard model with phonons it is necessary to 
consider two distinct representations of the superconductive Uq{su{2)) for different 
lattice sites. One (p^) is equal to p^, the other (p~) differs from p^ by a minus sign 
on the generators X^: 



pf(X+) = ±p,(X+), pf = ±p,(X-), pt{H) = p,{H). (4.4) 

For each lattice site / a sign a {I) G {1,-1} is chosen and the representation p+ or 
p~ is associated to it depending on whether a{l) = 1 or a{l) = —1 respectively. The 
local commutation relations are not affected by this choice. The sign will however 
be crucial for the global commutation relations. Hence, for the moment we will not 
specify a rule for assigning a representation to a given site, but we will see later ([4.29|) 



that sites corresponding to nearest neighbours must have opposite representations 
p"*" and p~. This is exactly what happens in the classical case 0. For orthogonal 
(square) lattices a choice of the sign which implements this condition is 

cr{l) = (-l)ll'l' (4.5) 

where ||/|| = Y.n=i is the length of the index I which labels the site /. 

For the moment we will choose an arbitrary ordering of the lattice sites. Choosing 
an ordering is necessary to be able to define a tensor product and hence to construct 
a global symmetry. According to the definition of the coproduct in Uq{su{2)) (see 



2.22| and footnote 4) 

A(X+) = e5"^ ®X+ + X+ ®e-^"*^, (4.6) 
A(X-) = e5"*-f^ O X- + X- ® e-^"^ = (A(X+))^ (4.7) 
A{H) = H^1 + 1^H, (4.8) 

where the deformation parameter is chosen to be g = e" and a is a complex para- 
meter to be determined by the commutation relations and through the representa- 
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tions pf we obtain the generators of global superconductive Uq{su{2)) 



KM 



,^^(0(A(A^-i)(X+)), 

,p:«(A(^-i)(x-)) = (;rW)t, 

,p^(0(A(^-i)(if)), 



(4.9) 
(4.10) 
(4.11) 



where is the number of lattice sites. Using ( [4.6H4.8|) these generators are computed 
to be 



I r<l r>l 

I r<l r>l 



and K^'^ = Y.^> 



(4.12) 
(4.13) 
(4.14) 



Local commutation relations The local part of the Hamiltonian commutes with 
the local generators 



if the following conditions hold 







(4.15) 



2A 



Muo^h^^^ = -. 

2 4 2 Mcj2 



(4.16) 
(4.17) 



Global commutation relations The fact that K^^^ commutes with 
given by ( p.23| ) is immediate. We must calculate 



(non—loc) 



(4.18) 



/ r<l r>l 



(2) 



{i<j) 'y 



^(Rijix^-Xj) piK-{pi-pj) t 



It can be seen that 



1 + 2Ki'\l - e-t) + 2n,|ni| ( cosh(-) - 1 



2sinh(^C^^ii ■ - S^)e-'^-^'+^^^^<^'-'^^\{A.19) 

(4.20) 
(4.21) 
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and, using equation ( |4.21| ), 



(4.22) 



+ (Sij{a]^ai^nji + rij-^a^^aii) - d^ia^^ai-^riii + nita]^ai|)) {e" + e " - 2). 
We introduce the abbreviation 

- - - (4.23) 



n e 



-a*K, 



Splitting the commutators, evaluating the expressions that are obtained by the use 
of ( [4 .191 - 1^^2^ ), and using the delta-functions which appear in ( |4.19| ), ( [4.2q ) and 
( [4.22|) to perform some of the sums, it can be seen that (|4. 18| ) becomes 



}{(+) Jj{non~loc) 



{i<i) 



Z,j ( 2 cosh ( -R,j ■ <l>C/ij - 2 cosh \^Rij ■ + -ct 



+ Zji (2 cosh {^Rij ■ *C^) - 2 cosh i^Rij ■ ^(^ + ^aj 



Zijie^"' - 1) + Zji (e-^^^-'^^^e-^" - 1 
%(e^--*^^e^"* - 1) + Z,,(e-^" - 1) 



I (i,j),i<l<j 



(4.24) 



There are 2 sums containing 6 fermionic operators, 4 sums containing 4 fermionic 
operators, and 2 sums containig 2 fermionic operators. These sums must all be 
separately 0, because they depend on different numbers of such operators and hence 
are linearly independent. Let's study the term containing aj|a]| — al-^Oji'- 



{i<3) 



(4.25) 
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The above sum can vanish only if each term with fixed i,j is separately 0, because 
there are no other terms which contain al|Ctj| "'^It'^ji- Therefore it is necessary that 
the expression between the square brackets is 0. For this reason we must require 

%e^-^'^-*^^e^"* + Zjie"5^'^*C?ig-^« = q (426) 
This is equivalent to the set of equations 

Zij = -Zji, (4.27) 

^lR,y$(;h^^a* ^ ^-^R^y$(;n^-la^ (4.28) 

which in turn imply — are nearest neighbours) 

a(t) = -a{j), (4.29) 

^-i{S-K)-pi-iK-pj+(;Rij-{xi-Xj) _ ^-i{^-K)-pj-iK.-pi+C,Rij-{xi-Xj) gg^ 

ne'^''''' n e-'^*^'^'= n e'^^'^'ne""*'''^ (4-31) 

r<i r>i,rj^j r<j,rj^i r>j 

^-^Ima^^^R^jScn^lRea _ ^-^R^y^(;h^-^Rea-^ ^ g_ (4.32) 

Eq. ( [4.29|) means that nearest neighbours must have opposite signs. As we have 
already anticipated, this means that in order for the global commutation relations 
to hold, it should be possible to see the lattice A on which the model is defined, as 
the sum of two lattices Ai, A2, such that nearest neighbours are always on different 
lattices. This gives a restriction on the possible lattices, e.g. a triangular lattice 
could not be chosen. 

Eq. ( [4.30|) implies k — $ = — k and hence 

2k = $. (4.33) 



This is one condition that must be satisfied for expression ( 4.25 ) to vanish. In 



particular it fixes the parameter of the transformation ( p.l7|) . It turns out, that the 
parameter has to be the same as the one used to transform the fermionic operators 
in the hamiltonian. 
Eq. (|;32D implies 

Rea = -Rij ■ $Ch. (4.34) 
This is the second condition which must be satisfied for expression ([4.25|) to vanish. 



It is important to notice that it is possible to fulfill this relation only if the ordering 
of the lattice sites is chosen to be the lexicographic one. So this imposes a first 
restriction on the ordering of the sites. 
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The strongest relation is ([4.31D — it depends crucially on the ordering chosen for 



the lattice sites. In order for ( |4.31| ) to hold it is necessary that 



n 

i<r<j 



(^) 



n e- 

i<r<j 



(4.35) 



Let's apply ( 4.21 ) to expand the exponentials. Then we obtain an expression of the 
type 

1 + 2(1 - e^) ^^'^ + . . . = 1 + 2(1 - ) + ■ ■ ■ (4-36) 

i<r<j i<r<j 

(Here the terms which are indicated with ". . ." are at least quadratic in the K^^^ 
and therefore are independent from the first-order terms which have been written.) 
Eq. ( 4.3(i| ) shows that in order for relation (4.35) to hold, it is necessary that 



e 2 



Re{a) = 0. 



But this would mean that the coproduct should be symmetric, and this is against 
the hypothesis that there is a true quantum symmetry. 

This shows that we must look for a condition on the ordering of the lattice sites, 
so that we do not need to require ( [4.35| ): There cannot be any site r which satisfies 
the condition i < r < j for any couple of nearest neighbours In other words it 
is necessary that if i,j is a couple of nearest neighbours then 



{i < r ^ j < r, i > r ^ j > r) Vr. 



(4.37) 



However, condition ( [4.37|) implies that the lattice A on which the Hubbard model 
is defined is one-dimensional, and that the 'normal' ordering of the sites is chosen, 
in which the sites are numbered from left to right in increasing or decreasing order. 

It can be verified immediately that condition ( [4.37|) is sufficient to guarantee 
that the sum with "X^z J2{i,j),i<i<j" is not present, because there is no longer any / 
satisfying i < I < j. In fact it can be verified with arguments similar to the ones 
used to study the term (|4.25|) that ( [4.37|) is also necessary for such a sum to vanish. 
Because of ( [4. 311 ) and ( [4.30| ) it is possible to combine the terms which contain the 
same products of fermionic operators. Further, it can be seen immediately, that 
conditions ( [4.29|) , ([4.33|) , (|4.34| ) and ([4.37|) are also necessary and sufficient for the 
sums with 4 and with 6 fermionic operators to vanish. Thus in order for the hamil- 
tonian (g^^ commute with the generators (^12| - gl4D of Uq{su{2)) 

the conditions ( [4.29|) , (|4737|) , (|4.3^^ ) , ( |4.34|) are necessary and sufficient. 
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5 Discussion 



The computation of the previous section has shown that the Hamiltonian H, given 
for a D-dimensional lattice in ( p. 21 - 3^2^) , commutes with the generators of (global) 



UgSM(2), q = e", provided that local conditions (^]T^, [4.17|) , global conditions (|4.33| , 



4.34| ) and either D = 1 or a = hold. (We will only consider Im(a) = here.) 



We are now in the position to comment on the symmetries^ of various hamiltonians 
derived from the Hubbard model: 

//nub ( |3.14[ ) coincides with ( |3.21| - pT^BD for the particular choice of parameter 



K = 0. The global conditions ([4.33| , ^4.34|) imply $ = 0, a = and hence q = 1. 
The local conditions ( [4.16| , \i.l7\) imply A = and /i = u/2. We conclude that the 



Hubbard model with phonons ( ^.141 ) has no true quantum symmetry — not even in 
the one-dimensional case; it has an ordinary SU(2) symmetry provided that there is 
no local electron-phonon coupling (A = 0). 

The Hamiltonian H ( |3.21| ) studied in the previous section and considered by 
P] is formally obtained from the Hubbard model with phonons by a Lang-Firsov 
transformation on the fermionic operators only. (See remark at the end of section ^) 
The essential difference between i^Hub and H is that -ffHub uses coordinates yi that 
commute with bj^ , hia while H is written in terms of new coordinates that commute 
with and (and not with hi„). To be able to compare the two Hamiltonians 
we have to relate the sets of coordinates. This is simply done by a Lang-Firsov 
transformation^ (see 3.15 - 3.17 ) 

Xi = U{R)mU~\ii). (5.38) 
The new coordinates are found to be 

Xi = yi + TiK^rii^, (5.39) 

i.e. the position of the ion at lattice site i is shifted according to the number of 
electrons at that site. Expressing H in terms of hi„ we find 



i i,a i \ / 



i \ {i<j) o- 

X (l + (e^^^- « - l)ni,_,) (l + (e-^^^-"^ - l)n,- + h.c) (5.40) 



^The Hamiltonian will commute with all elements of U<jSU(2) provided that its generators do 
so; this is equivalent to a full quantum symmetry (under quantum adjoint action), see section |^. 
^This observation is also supported by the choice of Ki^\Ki^\ 
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with a set of new parameters 

A' = X-Muj^M, (5.41) 

u = u - 2h\ ■ a + Muj\^ , (5.42) 

/i' = /i + /lA • K - 1 /2Muj'^Ti^K^. (5.43) 



The model given by ifq_sym ( ^-401 ) has a true quantum symmetry (in the one-dimen- 



sional case). The local conditions ( [4.16| , |4.17|) for quantum symmetry expressed in 
terms of the new parameters is 

A' = 0, /i' = u'/2. (5.44) 

There is apparently no local coupling to the phonons and the condition for symmetry 
is "half filling" as in the standard Hubbard model. 
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